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Abstract 

We study the Weyl representation of metaplectic operators associated 
to a symplectic matrix having no non-trivial fixed point, and justify a for¬ 
mula suggested in earlier work of Mehlig and Wilkinson. We give precise 
calculations of the associated Maslov-type indices; these indices intervene 
in a crucial way in Gutzwiller’s formula of semiclassical mechanics, and 
are simply related to an index defined by Conley and Zehnder. 
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1 Introduction 


In a remarkable paper M Mehlig and Wilkinson propose a simple derivation 
of Gutzwiller’s f2_ approximation 


PGuxJE) = Re J2 

po 
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for the oscillating part of the semiclassical level density for chaotic systems 
whose periodic orbits “po” are all isolated and non-degenerate ( T po is the prime 
period, v po an integer related to the Maslov index, A po the action, and S po 
the stability matrix). Mehlig and Wilkinson’s derivation heavily relies upon 
their observation that for any symplectic matrix S such that det(«S — ^ 0 

one has 
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where R U (S ) is the operator L 2 (M”) —» L 2 (M”) defined by 
/ 1 \ n MS) r . 

^*w=y jms-dJ '**™**™*^* 

where: 

• T(zo) is the Weyl Heisenberg operator associated to zo = (xo,Po) G IR 2n : 

T(z 0 )f(x) = e *(M-±<po,*o» /(x _ Xo) 

for any function / defined on M 2n ; 

• Ms is a symmetric matrix, associated to S by the formula 

M S = lJ(S + I)(S - I)- 1 (2) 

I being the 2nx2n identity matrix and J the standard symplectic matrix; 

• The integer v{S) corresponds to a choice of argdet(5 — I). 

In what follows we will write Mehlig-Wilkinson’s formula as a Bochner 
integral 

/ i \ n s) r 

MS) = (s) Vldet(5-/)| / «*<''—(3) 

The validity of Mehlig and Wilkinson’s derivation requires that -as these 
authors claim- R U (S) is one of the two metaplectic operators ±5 associated 
with the symplectic matrix S. To sustain this claim the authors refer to 
previous work of one of the authors E 5 they also claim that for all S, S' such 
that det(S — 7) / 0, det(S" — 7) / 0 and det (SS 1 — 7) ^ 0 their operators 
satisfy 

R u n(SS') = ±R v (S)R v ,{S). (4) 

The purpose of this paper is to fully justify Mehlig and Wilkinson’s state¬ 
ments. There are actually several options available. The a priori most direct 
strategy would be to use Howe’s theory m of operators with Gaussian ker¬ 
nels (the “oscillator semigroup” theory: see 0 for a review); another approach 
would be to use Hannabuss’ character theory [B] for contact transformations 
(also see the follow-up [I; to this paper for interesting applications to star- 
products). These two methods however both have, drawbacks. Howe’s theory 
needs quite a lot of technical prerequisites and would lead here to unneces¬ 
sarily long calculations; in addition it would not be very helpful for the study 
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of the sign ambiguity in © since this point is not really addressed in Howe’s 
work. On the other hand, Hannabuss’ machinery is quite abstract (it makes a 
heavy use of group character theory) and the use of this approach would per¬ 
haps have a tendency to obscure things. For these reasons we prefer a more 
straightforward line of attack, using standard Weyl calculus together with the 
theory of the metaplectic group as developed in [I]. This approach moreover 
has, as we will see, the overwhelming advantage of producing simple formulae 
relating the integer v in © to the usual Maslov index of the metaplectic group 
(we emphasize that v is not the Maslov index!). This is important, because 
Gutzwiller’s theory has been plagued since its very beginning by the question 
of how to calculate the indices v po appearing in the trace formula m , as wit¬ 
nessed by the abundant literature devoted to this delicate topic (see 17] where 
we discuss these issues and give a rather exhaustive list of references). 

Notations 

We denote by a the canonical symplectic form on the phase space = 

R” x R”: 

a(z,z') = (p,x r ) — (p,x) if z = ( x,p ), z' = (x'p') 
that is, in matrix form 

a(z, z') = ( Jz, z') , J = 

The real symplectic group Sp(n) consists of all linear automorphisms S of R^ n 
such that a(Sz,Sz') = cr(z, z') for all z,z'. It is a connected Lie group and 
iri(Sp(n)) is isomorphic to (Z,+). <S(R n ) is the Schwartz space of rapidly 
decreasing functions on R n and its dual S'( R n ) the space of tempered distri¬ 
butions. 

We will denote by Inert R the number of negative eigenvalues of a real 
symmetric matrix R. 

2 Prerequisites 

In this Section we briefly recall the main definitions and properties of the 
metaplectic group and of Weyl calculus we will need in the rest of this paper. 

2.1 Standard theory of Mp(n): Review 

The material of this first subsection is quite classical; see for instance im a ia 
and the references therein. 


0 / 
-I 0 
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Every S £ Mp(n ) is the product of two “quadratic Fourier transforms”, 
which are operators Sw,m defined on 5(M n ) by 


Sw,mf(x) = j e iW ( x,x '^ f (x')d n x' ( 5 ) 

where W is a quadratic form in the variables x, x' of the type 

W(x,x r ) = \{Px,x) — ( Lx,x') + \ (Qx' ,x') (6) 

with P = P T , Q = Q t , detL ^ 0. The integer m (“Maslov index”) appearing 
in © corresponds to a choice of argdetL: 

mi: = arg det L mod 2ir 

and to every W there thus corresponds two different choices of m modulo 4: 
if m is one choice, then m + 2 is the other, reflecting the fact that Mp(n) is 
a two-fold covering of Sp(n). The projection n : Mp(n ) —> Sp(n ) is entirely 
specified by the datum of each 7r(SV,m), and we have 7r(SV,m) = Sw where 

(x,p) = Sw{x',p') = d x W(x,x r ) and p' = — d x /W(x, x'). (7) 

Rewriting these conditions in terms of P,L,Q we get p = Px — L T x' and 
p' = Lx — Qx'\ solving these equations in x and p yields 

x = L~ 1 (p + Qx') and p = (PL~ l Q — L T )x' + PL~ x p> 

hence the projection S\y of S\y,m is just the free symplectic matrix 

o = r L~ l 1 

‘ bu PL-'Q - L t PL- 1 


generated by the quadratic form W. Note that if conversely S' is a free sym¬ 
plectic matrix 


S 


A B 
C D 


£ Sp(n ) 


det B 0 


(9) 


then S = Sw with P = DB ” 1 , L = R” 1 , Q = B _1 A. Observe that the free 
symplectic 2 n x 2 n matrices form a dense subset of S'p(n); this property will 
be used in the proof of Proposition II 01 

The inverse S^ ni = ( Sw,m )* of SV,m is the operator Sw*,m* where W*(x, x') 
—W(x',x) and m* = n — m, mod 4. 
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2.2 Heisenberg—Weyl operators 

The operators T(zq ) satisfy the metaplectic covariance formula: 

ST(z) = T(Sz)S (S = vr(5)) (10) 

for every S G Mp(n) and z. In fact, the metaplectic operators are the only 
unitary operators (up to a factor in S' 1 ) satisfying (1101) : 

For every S G Sp(n ) there exists a unitary transformation U in 
L 2 (M n ) satisfying HJ\) and U is uniquely determined apart from a 
constant factor of modulus one. 

The Heisenberg-Weyl operators moreover satisfy the relations 


T(z 0 )T(zi) = e ia ( z °’ zl ^T(zi)T(zo) 

(11) 

T{z 0 + Zl ) = e -^ Zl )f (zq)T ( Zl ) 

(12) 


as is easily seen from their definition. 

2.3 Weyl operators 

Let a w = a w (x,D ) be the Weyl operator with symbol a (which we always 
assume to belong to some suitable class, allowing the integrals to be viewed 
as distributions): 

a w f(x) = ) n fJ e i ^ x -^a(^x + y),p)f(y)d n yd n p- 
where / G S(M ra ); equivalently 

a w = (^r / a a (z 0 )T(z 0 )d 2n z 0 

where the “twisted symbol” a a is the symplectic Fourier transform F a a of a: 

a a (z) = F a a(z) = (£)»/ e~ icT ^ z ’ z '^ a(z')d 2n z'. 

The compose c w = a w o b w (when defined) is the Weyl operator with twisted 
Weyl symbol 

= {^) n ( a cr *a b a ) (13) 

where 

a b(z) = J e^ z ' u) a(z - u)b(u)d 2n u (14) 

(see for instance Littlejohn eh, Wong EH)- 
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2.4 Generalized Fresnel Formula 

We will use the following formula, generalizing the usual Fresnel integral to 
complex Gaussians. Let M be a real symmetric m x m matrix. If M is 
invertible then the Fourier transform of the exponential exp(i (Mx,x) /2) is 
given by the formula 

(^) m/2 J e- i{v ’ u) e^ Mu ' u) d m u= \detM\-^ 2 e^ ssnM e-^ M ^ v ) (15) 

where sgn M . the “signature” of M. is the number of > 0 eigenvalues of M 
minus the number of < 0 eigenvalues. 

For a proof of this formula see for instance .5], Appendix A. 


3 Discussion of the Mehlig—Wilkinson Formula 

In this Section we show that the Mehlig-Willkinson operators © indeed are 
metaplectic operators. We begin by giving two alternative expressions for 
these operators. 

3.1 Equivalent formulations 

We begin by remarking that the matrix 

M s = ^J(S + I)(S-I)- 1 

is symmetric; this immediately follows from the conditions 
S G Sp{n) <=> S T JS = J <=> SJS T = J. 

Notice that for every M with det(M — ^ J) ^ 0 the equation 

M= \J(S + I){S - I)- 1 
can be solved in S, yielding 

S = {M-\J)~ 1 {M + \J)- 

the relation S G Sp(n ) is then equivalent to M being real and symmetric. 
Lemma 1 Let S £ Sp(n) be such that det(S — I) ^ 0. The operator 

R U (S) = (—^ * [ e^ {MsZ ’ z) T(z)d 2n z (16) 

\2ttJ v /|det(5-/)|.y 
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can be written as 

R„(S)= (J^j *V|det(S-J)| J e-^ Sz ^T((S - I)z)d 2n z (17) 

that is, as 

MS) = zV|det(S-J)| j T(Sz)T(—z)d 2n z. (18) 

Proof. We have 

%J{S + I)(S - I)- 1 = + J{S - I)- 1 

hence, in view of the antisymmetry of J, 

{M s z, z) = (J(S - iMz, z) = a((S - M x z, z) 

Performing the change of variables z i—> (S — I)z we can rewrite the integral 
in the right-hand side of m as 

J e y MsZ ’ z) f(z)d 2n z = v/|det(5-/)| J e M’( s - I)z) T({S - I)z)d 2n z 

= vTdet {S-I)\ J e~i a ( Sz ’ z) T{(S - I)z)d 2n z 

hence m Taking into account the relation la we have 

T((S - I)z) = e^ Sz ’ z) f{Sz)T{-z) 
and formula m follows. ■ 

Remark 2 Formulae nw and Ol i OB* suggest that un could be a choice 
of ±argdet(S' — I). This is however not the case (see \2iA) in Proposition 0 
below); formula 12(A) will identify the integers with the Conley-Zehnder index. 

Corollary 3 We have R U (S ) = csSw,m with |cs| = 1. 

Proof. The operator R U (S) satisfies the metaplectic covariance relation 

MS)T(z) = T{Sz)MS ) 

as immediately follows from the alternative form CED of R U (S). On the other 
hand, a straightforward calculation using formula m shows that R U (S) is 
unitary, hence the claim. ■ 

Let us precise Corollary |3] by discussing the choice of the constant cs- 
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3.2 The case S — Si 


W,m 

We are going to show that the Mehlig-Wilkinson operators coincide with the 
metaplectic operators Sw,m when S = Sw and we will thereafter determine 
the correct choice for we will see that it is related by a simple formula to 
the usual Maslov index as defined in (3j. 

Let us first prove the following technical result: 

Lemma 4 Let Sw be a free symplectic matrix We have 

det(SV - I) = (-l) n det(.B) det (B~ l A + DB~ X - B~ l - (B T )~ 1 ) (19) 

that is, when S is written in the form JUb 

det(5vi/ - /) = (-l) n det(L _1 ) det(P + Q - L - L T ). (20) 

Proof. Since B is invertible we can write S — I as 


A-I 

B 


0 B 

C — (D — I)B~ 1 (A — I) O' 

C 

D-I 


I D-I 

B~ 1 (A — I) I 


and hence 

det(SV - I) = det(-P) det(C - (D - I)B~ l {A - I)). 

Since S is symplectic we have C — DB~ l A = —(B T )~ 1 (use for instance the 
fact that S T JS = SJS T = J) and hence 

C — (D — I)B~ 1 (A - /)) = B~ l A + DB" 1 - B~ l - ( B T )~ l ; 

the Lemma follows since det(— B) = (-T) n detb>. ■ 

Remark 5 The factorization m> shows in particular that ker(5 — I) is iso¬ 
morphic to ker (P + Q — L — L T ) (cf. I lty . Lemma 2.8, and proof of Lemma 
2.9). 

Let us denote by W xx the Hessian matrix of the function x i—> W(x,x), 
that is 

w xx = p + q-l-l t . 

We have: 

Proposition 6 Let S = Sw be a free symplectic matrix m and R U (S) the 
corresponding Mehlig-Wilkinson operator. We have R U (S ) = Sw,m provided 
that v is chosen so that 

v = m — Inert W xx mod 4 (22) 

in which case we have 

— argdet (S — I) =—u + n mod 2 
7r 
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Proof. Recall that we have shown that R U (S ) = csSw,m where cs is some 
complex constant with |cs| = 1. Let us determine that constant. Let 5 £ 
«S / (R”) be the Dirac distribution centered at x = 0; setting 

Cn> = (i) ~\/\ det(S — /)| 

we have, by definition of R U (S), writing zq = (xq,pq) in place of z = (x,p): 

R v (S)6(x) = C W)U J e^ MsZo ^e i ^ po ^-^ po ’ xo ^5(x-x 0 )d 2n z 0 
= C w ,u f e 3<^s(*j»o),(*jw)> e i<Po.*>$( x _ X0 )d 2n z 0 

hence, setting x = 0, 

Ru(S)8( 0) = C w ,u j e^ Ms( -°' po U 0 ’ po »6{-xo)d 2n z 0 
that is, since f S(—xo)d n xo = 1, 

RJS)S(O) = (—) , %V = [ q. (24) 

v ; v \27tJ y/\ det (5 — 1)\ J V ' 

Let us calculate the scalar product 

(M s (0,po), (0,p 0 )) = cr((5 - /) _1 0,p 0 ), (0,Po))- 

The relation (x,p) = (S — I) _1 (0,po) is equivalent to S(x,p ) = (x,p + po) that 
is to 

p + Po = x) and p = —d x >W(x,x). 

Using the explicit form m of W together with Lemma 0] these relations yield 

x = (P + Q-L- L T y l p, and p = (L - Q)(P + Q-L- L T y l p Q 


and hence 


{M s (0,po), (0,p 0 )) = - (W xx p 0 ,p 0 ) . 

Applying Fresnel’s formula CC33 we get 

J J(Ms(0, PO UO, PO )) d n po = ^y/\~fs g nW xx]detWxxl l/2. 

noting that 


1 

v / |det(S'w -I)| 


det L\ l / 2 \ det W xx \ 1 ^ 2 


(25) 
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(formula (1201) in Lemma 0} we thus have 

Ru(S w )5( 0) = (^) n/2 ^ e -^ sgnW -|detL| 1 / 2 . 

Now, by definition © of Sw,m j 

% m 5(0) = (^) n/2 i m - n/2 |detL| 1 /2 

hence i u e ~^ senWxx = i m_Tl / 2 . It follows that we have 
v — \ sgn W xx = m — mod 4 

which is the same thing as (@ since W xx has rank n. In view of m we have 
b arg det(5'u/ — /) = n + m + arg det W xx mod 2; 
formula (USD follows using (HI . ■ 

Let us digress for a while on the integers m and Inert appearing in 
formula ( 1221 ) and discuss them from the point of view of calculus of varia¬ 
tions. It is for this purpose useful to recall that in Gutzwiller’s formula (of 
which Mehlig and Wilkinson precisely want to give a new approach using the 
operators R U (S)) the symplectic matrix S is obtained from the monodromy 
matrix of an isolated Hamiltonian periodic orbit. Let us go a little bit further. 
Consider a Hamiltonian flow (f t determined by some time-dependent Hamil¬ 
tonian H = H(z,t ) defined on M 2n x R t , and let zq = (xo,Po) be such that 
4>t( z o) = ~o f° r some T > 0. The Jacobian matrices S(zo,t) = Dcj) t (zo) are 
symplectic and satisfy the “variational equation” 

^-S( Zo ,t) = JH"(zo,t)S(z 0 ,t) 
at 

where H” (zq, t) is the Hessian matrix D 2 H(cj) t (zo ), t) . When t varies from 0 to 
T the matrices S(zq, t ) describe a path in Sp(n) originating at the identity and 
ending at S(zq, T) (the “monodromy matrix”). Suppose that 5(^o, T) is a free 
symplectic matrix Sw ; then there exist p andp' such that (xo,p) = Sw(xo,p'), 
that is, expressing W in terms of P, L, Q as in © and using 0,P= ( P~L T )x 0 
and p' = (L — Q)x o, that is 

p-p' = (P + Q — L - L t )x 0 . 

It thus appears (see for instance mm) that Inert W xx is Morse ’sea order 
of concavity of the periodic orbit through ^o- 
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Remark 7 In m Piccione and his collaborators use the order of concavity, 
which they identify with Inert W xx , to investigate the Maslov and Morse indices 
for periodic geodesics. The topic is also discussed at some length in Muratore- 
Ginnaneschi m in connection with the study of Gutzwiller’s formulae using 
field-theoretical methods. 

Perhaps even more interesting is the relation between the index v and an 
index defined by Conley and Zehnder in J2J. Let us denote by Spo(n) the set 
of all S £ Sp(n) such that det(S — I) / 0. We have 

Spo(n) = Sp+(n ) U Sp-(n ) 

where S £ Sp±(n) if and only if ±det(S — I) > 0. The sets Sp±(n) are 
connected and every loop in Spo(n) is contractible in Sp(n). Consider now a 
path S : [0, T] \ —* Spo(n) going from the identity to S = S(T) ( S(T ) may 
be viewed, if one wants, as the monodromy matrix of a periodic Hamiltonian 
orbit corresponding to a time-dependent Hamiltonian). The index of Conley 
and Zehnder associates to the path S an integer p cz (S) only depending on 
the homotopy class (with fixed endpoints) of that path, and such that 

signdet(5 — I) = (—l)”“^oz(S) 

that is, equivalently, 

— argdet(S — I) = n — p-cz(S) mod 2. 

7r 

It follows from formula (USD in Proposition El that we have 

u = p cz {S) mod2 (26) 

and from formula (H2D in the same proposition that 

Pcz(S) = rn + n — Inert W xx mod 2. (27) 

3.3 The general case 

Recall that we established in Lemma El the equality 

det(SVi’ - I) = (- l) n det L~ x det(P + Q-L- L T ). (28) 

valid for all free matrices S\y £ Sp(n). Also recall that every S £ Mp(n) 
can be written (in infinitely many ways) as a product S = Sw.mSw,m' ■ We 
are going to show that Sy/,m and Sw',m' always can be chosen such that 
det (Syr — I) 7^ 0 and det(SV' — I) 0. For that purpose we need the 
following straightforward factorization result, which we nevertheless glorify by 
putting it into italics: 
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Lemma 8 Let W be given by B); then 

Sw,m = V-pML,mJV-Q (29) 

where 

V- P f{x) = e^ Px ’ x) f{x) , M Ltm f(x ) = z m y / | det L\f (Lx) , 
and J is the modified Fourier transform given by 

Jf( x ) = i ^) n/2 f e-^Kx'^x'. 

Proof. It is obvious using the explicit expression © of the quadratic form 
W (see |1|). ■ 

Let us now state and prove the first result of this section: 

Proposition 9 (i) Every S £ Mp(n) can be written as a product 

S = R V {S W )R V ,{S W ,). (30) 

(ii) The Mehlig-Wilkinson operators thus generate Mp(n). 

Proof. (ii) follows from (i) since the Sw,m generate Mp(n). To prove (i) 
let us write S = Sw,m.Sw',m’ and apply (El to each of the factors; letting 
P',L’, Q l define W l just as © is defined by P, L,Q we have 

S = V-pML,mJV_( P i + Q)ML' ,m' JV-Q'. (31) 

We claim that Sw,m and Sw',m' can be chosen in such a way that det (S\y—I) 7^ 
0 and det(5u/' — I) 0 that is, 

det(P + Q - L - L t ) 0 and det(P' + Q 1 - L 1 - L ,t ) ^ 0; 

this will prove the assertion in view of El. We first remark that the right 
hand-side of EB obviously does not change if we replace P' by P' + A and 
Q by Q — A where A € M. Choose now A such that it is not an eigenvalue of 
P + Q — L — L t and —A is not an eigenvalue of P' + Q' — L' — L' T then 

det(P + Q - XI - L - L t ) 0 and det(P' + A + Q' - L - L T ) 0 

and we have S = S\v 1 ,m. 1 Sw[,m' 1 with 

W\{x,x') = ^(Px,x) — (Lx,x') + \{{Q — X)x',x') 

W[(x,x') = \{{P' + A)x, x) — {L f x, x') + ^(Q'x', x'). 

■ 

So far, so good. But we haven’t told the whole story yet: there remains to 
prove that S £ Mp(n) can be written in the form R V {S) if det(5 — I) 0. 
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Proposition 10 Let S € Mp(n) be such that det(S' — I) 7 ^ 0. If S = 
R u {Sw)Rv'(S W ') then S = R U ^(S) with 

v(S) = v + v' + n — Inert(M + M') (32) 

the matrices M and M' being associated to S and S' by formula 11- 


Proof. A straightforward calculation using the composition formula and 
the Fresnel integral m shows that 


5 = 


1 

27T 


u+u'+^sgn(M+M') 


VI det(S w - I)(S W , - I)(M + M')\ 
where M and M 1 correspond to S\y and S\y by 0 and 

N = M — (Af + \J){M + -\J) 


e^ Nz ^T{z)d 2n z (33) 


We claim that 


det[(S w - I){S W f - I) (M + M’)\ = det(5 - I) (34) 

(hence M + M' is indeed invertible), and 

N= ±J(S + I)(S - I)- 1 = M S . (35) 

Formula (Era is easy to check by a direct calculation: by definition of M and 
M’ we have, since det J = 1, 


det[(S w ~ I)(S W , - I)(M + M’)\ = 

det[(5 w - /)(/ + (S w - I)- 1 + (. S w , - I)~ l ){S w , - I)} 


that is 

det[(Spv — I)(S\y — I)(M + M')\ = det(SVSV' — -0 

which is precisely (MI). Formula (E>51) is at first sight more cumbersome, and 
one might be tempted to use the oscillator semigroup calculations of Howe m 
at this stage. There is however an easier way out: assume that S = Sw",m"\ 
we know by Proposition El that we must have in this case 

N = ^J(SwSw + I){SwSw ~~ -0 1 


and this algebraic identity then holds for all S = SwSw' since the free sym- 
plectic matrices are dense in Sp{n). Formula (1331) can thus be rewritten 


S 


l \n jy +v ' + ±sgn(M+M') 

2^J VI det (S ~ W 


e U M sz,z)f( z )d 2n z 
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and formula El follows noting that if R is any real invertible 2 n x 2 n sym¬ 
metric matrix with q negative eigenvalues we have arg det R = qir mod 2vr and 
\ sgn R = n — q and hence 

-sgn(M + M') = n — Inert (M + M 1 ). 


4 Concluding Remarks 

We have justified Mehlig and Wilkinson’s claim that the metaplectic operators 
corresponding to symplectic matrices with no eigenvalues equal to one can be 
written in the form 0 ; we have in addition shown that every metaplectic 
operator can be written as the product of exactly two such operators. There 
are however still interesting open problems. It would be interesting to relate 
the index v appearing in ED to the cohomological Maslov index on Mp(n ) we 
constructed in HH: this would certainly lead to simpler -or at least more 
tractable- calculations for the indices intervening in Gutzwiller’s formula, as 
already demonstrated in our previous paper 0 where we examined the Maslov 
index of the monodromy matrix associated to a periodic Hamiltonian orbit. 
As shown by El a related mathematical problem would be to express the 
Conley-Zehnder index in terms of Leray’s index studied in de Gosson |2 d . 

We hope to come back to these important and interesting questions in a 
near future, together with applications to various trace formulae. 
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